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摘 要 : 为 了 快速 有 效 检测 南 疆 地 区 典型 土壤 (沙壤土 ) 的 盐分 含量 变化 ,利用 光谱 仪 和 电导 仪 测 得 
南 性 阿拉 尔 市 红 囊 种 植 区 盐 渍 土 近 红 外 高 光谱 和 电导 率 数 据 ,基于 7 种 不 同 光 谱 预 处 理 方法 和 2 
种 特征 波长 选择 算法 ,分 别 建立 多 元 线性 回归 (MLR) 和 偏 最 小 二 乘 回归 (PLSR) 的 土壤 盐分 监测 
模型 。 结 果 表 明 :7 种 预 处 理 方法 中 , 归 一 化 ,多 元 散射 ,变量 标准 化 和 一 阶 导 数 能 够 有 效 提高 土壤 


盐分 的 预测 模型 精度 。 基 于 多 元 逐步 回归 (SMR ) 


波长 选择 方法 的 多 元 线性 回归 (SMLR ) 模型 的 


R, <0.948 9,RPD <6. 294 9,RMSEP >0.435 6; 基 于 连续 投影 算法 (SPA) 的 多 元 线性 回归 (SPA- 
MLR) # #1 8) R? <0.956 8,RPD <6. 922 1,RMSEP >0.361 6, 预 测 结果 要 优 于 偏 最 小 二 乘 回归 


(PLSR ) 模 型 ,其 中 基于 归 


化 处 理 后 的 SMLR 和 SPA-MLR WW HW ER AER ARIA RS = 


0.979 2,RPD =9. 907 8 , RMSEP =0. 287 6 和 R,, =0. 980 5,RPD =10. 50,RMSEP =0. 278 3 ,而 且 
第 选 的 特征 波长 较 少 。 说 明 归 一 化 是 更 有 效 的 光谱 预 处 理 方 法 ,多 元 线性 回归 ( MLR) 更 适合 建立 


南 疆 典型 沙壤土 盐分 食量 的 预测 模型 。 
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土壤 盐 渍 化 是 新 疆 南 部 地 区 土地 退化 的 主要 问 
题 ,对 南 吐 农业 可 持续 发 展 以 及 生态 环境 的 保护 至 
关 重 要 "| 。 传 统 土壤 盐 渍 化 监测 ,依赖 于 土壤 理化 
分 析 方 法 ,不 能 快速 有 效 的 反映 区 域 土壤 盐 省 化 动 
态 信息 “ ,近年 来 通过 土壤 光谱 信息 估算 土壤 中 全 
氮 ”、 水 分 “、 有 机 质 等 研究 较 多 ,利用 高 光谱 预 
测 土壤 盐 渍 化 定量 分 析 也 逐步 展开 , 双 感 技术 利用 
土壤 的 盐分 光谱 特征 ,能够 快速 .实时 的 获得 盐 渍 化 
土壤 的 盐分 信息 “ … ,因此 ,研究 土壤 盐分 的 光谱 特 
征 ,建立 基于 光谱 分 析 的 土壤 盐分 定量 反 演 模型 是 
一 项 重要 的 工作 。 

土壤 盐 渍 化 的 高 光谱 定量 反 演 模型 是 土壤 盐 汗 
化 高 光谱 遥感 监测 的 重要 研究 内 容 之 一 。 目 前 ,已 
有 学 者 利用 不 同 的 研究 方法 对 不 同 地 区 的 盐 渍 化 土 
壤 盐 分 进行 了 高 光谱 定量 反 演 。 常 用 的 建 模 方法 包 
括 偏 最 小 二 乘 回归 、 多 元 线性 回归 和 主 成 分 回归 等 。 
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GOLDSHLEGE NI 等 通过 获取 番 若 地 高 光谱 反射 
率 与 土壤 盐分 含量 ,使 用 偏 最 小 二 乘法 建立 预测 模 
型 ,取得 了 较 好 的 预测 精度 , 其 中 预测 均 方 根 误差 
RMSE 是 0.6% ,决定 系数 尼 Æ 0.87, NETO OR”! 
等 在 巴西 半 干 旱地 区 使 用 土壤 电导 率 进 行 土壤 盐 渍 
化 评价 ,对 比 使 用 主 成 分 回归 ,多 元 线性 回归 和 偏 最 
小 二 乘法 建立 预测 模型 ,其 中 偏 最 小 二 乘法 获得 了 
较 好 的 预测 精度 , RMSE 是 0. 22,RPD 是 2.21, % 
杰 等 "中 通过 分 析 南 驻地 区 土壤 的 高 光谱 数据 与 盐 
分 含量 关系 ,基于 连续 统 去 除法 所 建 模型 的 相对 分 
析 误 差 (RPD) 达 2.5 WE. Jaa SE" 建立 了 内 
蒙古 河套 灌区 土壤 盐分 与 高 光谱 数据 的 偏 最 小 二 乘 
回归 (PLSR ) 模型 ， 验 证 数据 的 决定 系数 ( 尼 ) 为 
0.73 ,预测 方差 比 为 1.79。 丁 建 丽 等 通过 分 析 新 
疆 塔 里 木 盆 地 北 缘 渭 干 河 一 库 车 河 三 角 洲 绿洲 盐 汪 
化 土壤 光谱 信息 与 盐分 关系 ,发 现 利 用 综合 光谱 指 


国家 自然 科学 基金 项 目 (61501314 ,41561088 ) ;塔里木 大 学 校长 基金 项 目 (TDZKQN201614) 
j 安 人 ,塔里木 大 学 讲师 ,研究生 ,主要 从 事 遥 感 与 数字 农业 技术 研究 .Email:wujiang0156@ 163. com 
铁 成 (1982 - ) 男 ,黑龙 江 绥化 人 ,塔里木 大 学 副教授 ,比利时 列 日 大 学 博士 研究 生 ,了 


,主要 从 事 遥 感 与 数字 农业 技术 研究 . 


干旱 区 地 理 


数 来 监测 土壤 盐 渍 化 是 非常 有 效 的 ,预测 的 RMSE 
为 0.009。 段 鹏 程 等 扩 通过 对 新 疆 阜 康 500 水 库 下 
游 的 盐 渍 化 土壤 实地 定点 取样 和 光谱 测量 ,建立 了 
土壤 反射 光谱 与 盐分 含量 之 间 的 多 元 线性 回归 预测 
模型 ,在 不 同类 型 的 土壤 含 盐 量 预测 模型 中 ,一 阶 
导数 建立 的 模型 预测 效果 最 优 ,R 均 超过 0. 983, 
李 亚 莉 等 '… 以 新 疆 典 型 盐 渍 化 灌区 为 研究 对 象 , 获 
取 不 同 质地 的 土壤 水 盐 含量 光谱 曲线 ,采用 偏 最 小 
二 乘 回 归 方 法 构建 水 盐 反 演 模型 ,模型 的 决定 系数 
R 均 大 于 0.65。 

尽管 已 有 研究 利用 土壤 电导 率 和 含 盐 量 指标 ， 
取得 了 相对 理想 的 反 演 效果 , 但 不 同 区 域 的 土壤 结 
构 .成 分 差异 显著 ,不 同类 型 土壤 的 高 光谱 响应 和 敏 
感 波段 也 具有 一 定 的 差异 性 ,不 同方 法 建 模 的 预测 
效果 也 有 很 大 不 同 。 因 此 ,本 文 在 总 结 前 人 研究 基 
础 上 ,以 新 疆 阿 拉 尔 沙壤土 为 研究 区 ,依据 高 光谱 数 
据 与 土壤 电导 率 实 测 数据 ,对 获取 的 土壤 光谱 数据 
进行 多 元 散射 校正 矢量 归 一 化 ` 变 量 标准 化 、 移 动 
SER SG 平滑 一 阶 导 和 二 阶 导 等 不 同 预 处 理 ,分 别 
使 用 多 元 逐步 回归 (SMR ) 与 连续 投影 算法 (SPA ) Sif 
选 土壤 电导 率 敏 感 特征 波长 ,对 比分 析 多 元 线性 回 
归 (MLR) 和 偏 最 小 二 乘 回归 法 (PLSR ) 建立 土壤 电 
导 率 检测 模型 的 精度 和 效率 ,由 在 筛选 和 建立 一 种 
适合 该 区 域 的 土壤 盐分 预测 方法 和 模型 ,为 利用 遥 
感 技术 快速 .无 损 监测 新 疆 南 部 地 区 土壤 盐 渍 化 提 
供 理论 依据 和 实用 方法 。 


1 材料 与 方法 


1.1 数据 采集 

试验 区 位 于 新 疆 阿 拉 尔 市 ,81°13’3”E, 40°34" 
45”N, 束 园 土壤 类 型 为 沙壤土 。 分 别 采集 土壤 表层 
盐分 已 结晶 的 表层 土 和 0 ~20 cm 的 正常 土壤 , 按 
1:99,2:98 ,…,99:1 的 比例 混合 共 获得 142 份 不 同 
盐分 的 土壤 样本 , 带 回 实验 室 自然 风干 , 磨 细 过 筛 ， 
用 于 光谱 及 理化 分 析 测 试用 。 因 为 土壤 电导 率 已 成 
为 国际 上 广泛 使 用 的 盐分 表示 方式 | ,本 研究 使 用 
电导 率 代表 土壤 盐分 含量 状况 ,电导 率 采 用 水 土 比 
为 5:1 ,利用 DDS 一 307 电导 率 仪 直接 进行 测定 。 光 
谱 采 用 Zolix Gaia Sorter 近 红 外 成 像 高 光谱 仪 , 光 谱 
测定 范围 900 ~ 1 700 nm( 实 际 测量 到 1 750 nm), 
光谱 分 辩 率 5 nm ,光谱 采样 点 4 nm, 共 256 个 波段 
测定 获取 高 光谱 数据 ,每 次 测量 前 进行 标准 白板 校 
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正 ,每 个 样 点 重复 采集 10 条 光谱 曲线 , 取 平 均值 作 
为 该 样本 的 光谱 反射 率 ,剔除 了 2 份 异常 样本 后 
(光谱 曲线 缺失 ) , 共 获 得 140 个 有 效 样本 ,如 表 1 
是 根据 SPXY 算法 "挑选 出 的 100 份 样品 作为 定 
标 集 ,其 他 40 份 土壤 样品 作为 预测 集 。 
光谱 曲线 通过 异常 剔除 和 算数 平均 值 处 理 后 ， 
获取 140 个 样本 的 光谱 吸收 谱 ( 图 1)。 
1.2 光谱 数据 预 处 理 与 建 模 方法 
1.2.1 光谱 数据 预 处 理 方法 与 特征 波长 选择 方法 
由 于 近 红 外 光谱 区 的 吸收 强度 弱 ,光谱 信 噪 比 低 ， 
谱 带 易 重 县 ,为 了 消除 仪器 噪声 等 对 光谱 数据 的 影 
响 以 及 准确 分 析 土 壤 样 品 中 的 含 盐 量 ,有 必要 对 土 
壤 原 始 光 谱 进 行 预 处 理 以 提高 检测 的 精度 。 对 原始 
光谱 数据 分 别 进行 了 多 元 散射 校正 (MSC) .矢量 归 
一 化 (VN) 变量 标准 化 (SNV) 移动 平滑 (MA) SG 
平滑 (SG) 一 阶 导 (1 - Der) 和 二 阶 导 (2 - Der) 等 预 
处 理 , 其 中 MSC,SNV,VN 可 以 校正 样品 因 颗 粒 散射 
而 引起 的 光谱 误差 ,VN 可 以 较 好 的 改善 样品 稀 玻 


时 的 预测 结果 ,平滑 处 理 可 以 有 效 平滑 高 频 噪声 , 提 
高 信 噪 比 ,一 阶 导 数 和 二 阶 导数 分 别 用 于 消除 光谱 


中 基线 的 平移 和 漂移 ,可 有 效 消除 其 他 背景 干扰 ,分 


表 1 土壤 校正 集 和 预测 集 盐分 含量 统计 
Tab.1 Statistics of soil electrical conductivity in 


calibration and validation set 


数据 集 ”样本 数 ”最 小 值 RAE 平均 值 ”标准 差 

校正 集 100 0.2720 10.5400 5.6561 2.8131 

预测 集 40 1.1440 9.8700 1.1440 2.5422 
2 E 


吸光 度 
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图 1 土壤 原始 光谱 曲线 图 


Fig.1 Original spectrum of soil 
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土壤 电导 率 敏 感 特征 波长 选择 分 别 采用 多 元 逐 
步 回归 (SMR ) 与 连续 投影 算法 (SPA) 算法 ,其 中 
SMR 算法 可 集成 在 MLR 建 模 算 法 中 ,构建 多 元 线 
性 回归 模型 (SMLR ) 。 
1.2.2 建 模 方法 及 评价 指标 SMLR 是 在 多 元 线 
性 回归 的 基础 上 发 展 起 来 的 ,逐步 回归 按照 全 
局 自 变量 对 因 变 量 的 作用 、 显 著 程 度 或 贡献 率 ,逐个 
引入 到 回归 方程 中 ,剔除 对 因 变 量 不 显著 的 自 变 
ROO) ,建立 多 元 线性 回归 模型 ,对 因 变 量 进行 预 
测 。 

PLSR 吸取 了 主 成 分 分 析 中 提取 主 成 分 的 思想 ， 
能 够 简化 数据 结构 , 主 成 分 之 间 不 但 相关 性 最 小 ,而 
且 与 因 变 量 的 相关 性 最 大 ,同时 ,能够 克服 主 成 分 分 
析 对 自 变 量 有 较 强 解释 的 能 力 , 有 效 提取 对 系统 解 
释 性 最 强 的 综合 变量 ,提高 模型 的 估 测 能 力 。 因 而 ， 
本 研究 分 别 采用 多 元 线性 回归 和 偏 最 小 二 乘 回 归 方 
phe? 2 构建 土壤 电导 率 反 演 模型 ,预测 光谱 特征 值 
和 土壤 电导 率 之 间 的 关系 。 

预测 模型 采用 外 部 验证 ,对 其 稳定 性 和 预测 精 
度 进 行 评价 ,采用 建 模 决定 系数 (Ri? ) 、 建 模 均 方 根 
误差 (RMSEE) ,预测 决定 系数 (RR?) .预测 均 方 根 
误差 (RMSEP) .相对 分 析 误 差 (RPD ) 为 精度 评价 指 
标 ,R ”和 RPD BK, RMSEP 越 小 ,说 明 预 测 效果 
越 好 。 


2 数据 处 理 结果 与 分 析 


2.1 土壤 电导 率 多 元 线性 回归 模型 
对 原始 光谱 分 别 进行 VN 、MSC、SNV SG 等 预 
处 理 , 在 SMR 特征 波长 选择 的 基础 上 ,分 别 建立 多 
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元 逐步 线性 回归 (SMLR ) 模 型 ,并 进行 验证 。 首 先 
将 样品 电导 率 作为 因 变 量 , 设 定 模型 入 选 变量 的 下 
检验 概率 为 0. OS ,模型 剔除 变量 的 下 检验 概率 为 
0.1 ,将 其 对 应 的 土壤 原始 光谱 及 其 预 处 理 数据 作为 
自 变量 ,建立 SMLR 模型 ,其 预测 结果 如 表 2。 基 于 
7 种 预 处 理 方 法 建 模 方 程 的 Roa 和 预测 R,, 都 很 
高 ,而 且 都 大 于 交叉 验证 的 RR 。 其 中 , 归 一 化 ,多 
元 散射 校正 ,变量 标准 化 的 预测 精度 相对 较 好 且 优 
于 原始 光谱 ,预测 的 RMSEP 在 0.287 6 Ail 0. 331 2 
之 间 ,R,, 在 0.969 0 和 0.979 2 之 间 , RPD 在 
8. 104 4 和 9.907 8 之 间 , 其 中 归 一 化 处 理 后 的 预测 
精度 最 高 (RMSEP =0. 287 6,R ”=0.979 2,RPD = 
9.907 8). 

在 原始 光谱 预 处 理 基础 上 ,使 用 SPA 筛选 的 特 
征 波 长 作为 自 变量 ,建立 SPA-MLR 多 元 线性 回归 模 
型 ,对 应 的 土壤 电导 率 作为 因 变 量 , 预测 结果 如 表 
3。 对 于 7 种 预 处 理 方法 , 建 模 方程 的 Rs ”和 预测 
R 都 很 高 ,而且 都 大 于 交叉 验证 的 R, 。 移 动 平 
滑 和 二 阶 导 数 的 预测 精度 较 差 ,其 它 预 处 理 的 预测 
精度 相对 较 好 且 优 于 原始 光谱 ,模型 的 RMSEP 在 
0. 267 5 和 0.3079 之 间 ,R,, 在 0.968 9 和 0.9805 
之 间 ,RPD 在 8.1305 和 11.178 之 间 , 其 中 一 阶 导 
数 预测 精度 最 高 (RMSEP =0. 267 5 ,R,,” =0. 980 5, 
RPD =11.178)。 尽 管 经 过 MSC,SNV,1 - Der 和 2 - 
Der 预 处 理 后 SPA 筛选 的 变量 较 多 ,但 和 SMLR 的 
预测 结果 比较 ,SPA-MLR 预测 的 RMSEP 降低 ,Rs 
Al RPD 提高 ,说 明 SPA-MLR 的 预测 效果 略 好 于 
SMLR 方法 。 

2.2 土壤 电导 率 偏 最 小 二 乘 回归 模型 
原始 光谱 预 处 理 基 础 上 ,直接 使 用 PLSR 建 立 


表 2 基于 逐步 回归 SMR 的 土壤 电导 率 多 元 线性 回归 模型 (SMLR ) 
Tab.2 Soil electrical conductivity multiple linear regression model based on SMLR 
a 特征 波 定 标 交叉 验证 预测 
数据 变换 - 5 5 
段 个 数 RMSEE Rea RMSECV Ra RMSEP Ra RPD 
原始 光谱 4 0.477 4 0.941 0 0.505 0 0.934 2 0.388 6 0.953 2 6.579 8 
归 一 化 3 0.423 9 0.948 4 0.442 0 0.943 9 0.287 6 0.979 2 9.907 8 
移动 平滑 5 0.380 5 0.964 2 0.406 3 0.959 2 0.362 0 0.956 2 6.805 1 
多 元 散射 校正 3 0.415 5 0.952 8 0.435 0 0.948 3 0.331 2 0.969 0 8. 104 4 
变量 标准 化 3 0.397 0 0.956 9 0.416 6 0.952 6 0.331 4 0. 969 4 8.1549 
卷 积 平滑 4 0.424 6 0.955 0 0.446 6 0. 950 3 0.388 8 0.9515 6.462 0 
一 阶 导 数 5 0.474 2 0.937 9 0.506 4 0.929 3 0.425 6 0. 948 9 6.294 9 
二 阶 导数 6 0.419 2 0.950 6 0.458 9 0.941 0 0.435 6 0.949 4 6.325 7 


干 早 


土壤 电导 率 回 归 模 型 ,预测 效果 如 表 4。 其 中 主 成 
分 数 表示 从 土壤 光谱 曲线 中 提取 的 用 于 解释 土壤 电 
导 率 的 成 分 个 数 ,对 于 建 模 所 需 提取 的 主 成 分 个 数 ， 
可 以 通过 交叉 有 效 性 检验 来 确定 。 基 于 7 种 预 处 理 
建 模 的 R,, ”和 预测 R, 都 很 高 ,而 且 都 大 于 交叉 验 
证 的 R,。 除 了 卷 积 平滑 的 预测 精度 相对 较 差 ,其 
他 处 理 后 的 模型 预测 值 都 优 于 原始 光谱 ,预测 的 
R a 在 0.9543 和 0.9737 之 间 ,RMSEP 在 0.322 8 
和 0.400 4 之 间 ,RPD 在 6.679 8 和 8.838 1 之 间 ， 
说 明 在 光谱 不 同 预 处 理 的 模型 预测 精度 存在 较 大 差 


异 。 从 预测 的 Ra 来 看 归 一 化 ,移动 平滑 ,多 元 散 
射 校正 ,变量 标准 交 化 ,一 阶 导 数 ,二 阶 导 数 都 大 于 
0. 95 ,预测 效果 均 较 好 ,其 中 归 一 化 处 理 后 的 精度 最 
高 (RMSEP = 0. 322 8,R,, = 0. 973 7, RPD = 
8.838 1), 

原始 光谱 预 处 理 后 ,使 用 SPA 算法 筛选 敏感 波 
长 建立 偏 最 小 二 乘法 回归 模型 (SPA-PLSR) ,SPA Sif 
选 的 特征 波长 作为 自 变 量 , 对 应 的 土壤 电导 率 作 为 
因 变 量 ,预测 结果 如 表 5。 对 卷 积 平 滑 处 理 后 的 光 
谱 建 立 SPA-PLSR 模 型 ,没有 通过 0. 05 的 水 平 检 


val 


R3 基于 SPA 选择 特征 波长 的 土壤 电导 率 多 元 线性 回归 模型 (SPA-MLR ) 
Tab.3 Soil electrical conductivity multiple linear regression model based on SPA (SPA-MLR) 


P 特征 波 定 标 交叉 验证 预测 
段 个 数 RMSEE Rea RMSECV Ra RMSEP Rar RPD 
原始 光谱 5 0. 419 6 0. 954 6 0. 447 9 0. 948 2 0.311 8 0. 970 7 8. 444 5 
归 一 化 7 0.4147 0. 950 7 0.449 7 0. 942 1 0. 278 3 0. 980 5 10. 501 0 
移动 平滑 6 0. 392 7 0.961 8 0.424 7 0. 955 4 0.361 6 0. 956 8 6.922 1 
多 元 散射 校正 14 0.375 0 0. 961 6 0. 437 3 0. 948 1 0. 286 3 0.977 0 9.791 1 
变量 标准 化 14 0.372 2 0. 962 2 0.430 6 0. 949 5 0.304 9 0.974 1 9, 152 4 
卷 积 平滑 5 0. 405 0 0.959 1 0. 431 3 0. 953 7 0. 307 9 0. 968 9 8. 130 5 
一 阶 导数 18 0. 366 7 0. 963 1 0. 453 4 0. 944 0 0.267 5 0. 980 5 11.178 0 
二 阶 导数 12 0. 399 3 0. 955 2 0. 465 8 0. 939 5 0.341 6 0. 969 6 8. 227 3 
表 4 土壤 电导 率 偏 最 小 二 乘 回 归 模 型 (PLSR ) 
Tab.4 Soil electrical conductivity partial least squares regression model ( PLSR ) 
pan REY 难 证 iy ae pe sy 
数据 变换 定 标 交叉 验证 - 预测 主 成 份 
RMSEE Ri RMSECV Ra RMSEP Ra RPD 

原始 光谱 0. 6917 0.874 1 0.702 0 0. 871 1 0.556 9 0.916 0 5. 087 9 4 

归 一 化 0. 397 1 0.954 8 0. 448 7 0.941 9 0. 322 8 0. 973 7 8. 838 1 6 
移动 平滑 0.383 0 0.963 7 0.434 8 0. 953 4 0.365 7 0.954 3 6. 679 8 7 

多 元 散射 校正 0. 408 7 0. 954 3 0. 437 0 0. 947 8 0. 344 7 0. 966 7 7.852 6 4 

变量 标准 化 0. 406 5 0. 954 8 0. 434 5 0. 948 4 0. 343 5 0. 967 0 7.8775 4 

卷 积 平滑 0. 673 5 0. 884 9 0.7129 0. 872 1 0.6157 0. 884 5 4.489 4 4 
一 阶 导 数 0.356 8 0.965 1 0.491 3 0. 934 9 0. 337 8 0. 969 3 8. 850 3 7 
二 阶 导 数 0.386 3 0.958 1 0.450 0 0.942 2 0. 400 4 0.958 7 7. 176 5 5 

RS 基于 SPA 选择 特征 波长 的 土壤 电导 率 偏 最 小 二 乘 回归 模型 (SPA-PLSR ) 
Tab.5 Soil electrical conductivity partial least squares regression model based on SPA (SPA-PLSR ) 
a 特征 波 定 标 交叉 验证 预测 主 成 份 
段 个 数 RMSEE Rar RMSECV Ra” RMSEP Ra RPD 

原始 光谱 5 0.419 6 0.954 6 0.4479 0.948 2 0.311 8 0.9707 8.444 5 5 
归 一 化 7 0.577 6 0.903 2 0.601 2 0.895 3 0.449 2 0.950 5 6.832 8 3 
移动 平滑 6 1.095 5 0.684 9 1.145 8 0.661 4 1.051 5 0.589 4 2.375 1 3 
多 元 散射 校正 14 0.412 0 0.953 6 0.463 3 0.941 4 0.355 6 0.963 9 7.793 3 5 
变量 标准 化 14 0.422 1 0.951 3 0.448 7 0.945 0 0.370 9 0.960 8 7.206 0 4 
一 阶 导 数 18 0.424 9 0.950 3 0.454 4 0.943 2 0.324 2 0.970 1 8.563 1 4 
二 阶 导 数 12 0.472 1 0.937 1 0. 494 6 0.930 7 0.446 8 0.947 7 6. 309 3 4 


ChinaXiv 合 作 期 刊 


王涛 等 :MLR 和 PLSR 的 沙壤土 盐分 含量 光谱 检测 对 比 研究 


验 , 因 此 没有 建立 预测 模型 。 由 表 5 可 以 看 出 ,对 于 
6 种 预 处 理 方法 ,经 过 SPA 特征 波长 选择 后 的 建 模 
预测 精度 都 略 低 于 原始 光谱 ,尽管 多 元 散射 校正 和 
一 阶 导 数 的 预测 精度 相对 较 好 ,但 是 所 选择 的 特征 
波长 较 多 ,计算 量 较 大 。 这 可 能 因为 经 过 SPA 选择 
后 人 选 PLSR 的 主 成 分 与 直接 使 用 PLSR 的 主 成 分 
相 比 较 数 目 明显 减少 ,丢失 了 部 分 有 效 信息 导致 模 
型 的 预测 精度 有 所 下 降 。 

比较 表 4 和 表 5 发 现在 光谱 预 处 理 前 提 下 进行 
SPA 特征 波长 选择 后 建立 的 PLSR 模型 的 预测 精度 
要 低 于 直接 使 用 PLSR 模型 ,但 是 基于 原始 光谱 的 
SPA-PLSR 模型 的 预测 精度 高 于 PLSR 模型 。 分 析 
原因 经 过 SPA 选择 后 ,通过 选择 特征 波长 来 代表 全 
波段 进行 建 模 本 身 就 会 丢失 一 些 有 效 信息 ,在 通过 
PLSR 建 模 筛 选 主 成 分 后 ,有 效 信息 又 会 丢失 一 些 ， 
导致 模型 的 预测 精度 降低 。 比 较 表 3 和 表 5 发 现在 
光谱 预 处 理 前 提 下 建立 的 SPA-MLR 模型 精度 高 于 
PLSR ,但 基于 原始 光谱 建立 的 模型 预测 精度 却 是 一 
样 的 。 分 析 原 因 MLR 不 会 筛选 主 成 分 ,有 效 信息 相 
对 PLSR 丢失 较 少 ,因此 精度 相对 较 高 。 因 为 原始 
光谱 存在 噪声 导致 经 过 SPA 筛选 特征 波长 后 与 
PLSR 出 现 得 选 一 致 的 情况 ,2 种 模型 预测 精度 一 
致 ,从 而 也 说 明 ,光谱 在 没有 经 过 预 处 理 进 行 建 模 
时 ,因为 受到 噪声 的 影响 , 建 模 算法 失效 ,预测 结 
没有 了 可 比 性 。 

综 上 ,基于 归 一 化 处 理 后 使 用 多 元 逐步 回归 
(SMLR) 和 基于 连续 投影 算法 的 多 元 线性 回归 模型 


1 人 二 X 
SMLR 预测 的 散 点 图 
8 上 
& 6f 
= 
RMSEP=0.287 55 
R?=0.979 21 
2 p RPD=9.907 8 
0 1 , f 1 1 
0 2 4 6 8 10 
实测 值 /% 


K2 基于 归 一 化 预 处 到 


(SPA-MLR ) 的 预测 精度 相对 较 好 ,并 且 筛 选 的 特征 
波长 也 较 少 ,如 图 2 所 示 。 


3 讨论 


不 同 算法 的 预测 结果 比较 如 表 6 所 示 ,不同 变 
换 形式 对 土壤 光谱 信息 有 一 定 影响 ,进而 影响 佑 测 
模型 的 精度 。 

本 文 基于 不 同 光 谱 预 处 理 ,采用 SPA 选择 回归 
变量 建立 的 MLR 模型 的 预测 精度 要 优 于 SMR ,其 
原因 可 能 是 SPA 在 选择 特征 波长 的 能 力 上 要 比 
SMR 强 , 且 SPA 选择 的 特征 波长 要 多 于 SMR ,多 个 
相关 波长 能 够 较 好 的 反应 土壤 有 机 质 的 光谱 信息 ， 
从 而 可 以 较 好 估算 土壤 有 机 质 含量 ,但 是 SPA 选择 
的 回归 变量 较 多 ,计算 比较 耗 时 ,SMR 建 模 仅 有 少 
量 波 长 的 光谱 能 进入 模型 ,而 且 这 些 波 长 之 间 仍 然 
存在 一 定 的 共 线 性 ,模型 对 于 原始 光谱 信息 利用 较 


少 。 


直接 使 用 PLSR 方法 建 模 的 预测 精度 要 优 于 
SPA-PLSR 模型 。 直 接 使 用 PLSR ,可 入 选 的 主 成 分 
个 数 要 多 于 SPA-PLSR ,包含 了 更 多 的 有 效 信息 ,分 
析 原 因 是 PLSR 的 基本 作法 是 首先 在 自 变量 集中 提 
出 第 一 成 分 ,然后 建立 因 变 量 与 第 一 成 分 的 回归 方 
程 ,如 果 回 归 方 程 已 达到 满意 的 精度 , 则 算法 中 止 。 
否则 继续 第 二 对 成 分 的 提取 ,直到 能 达到 满意 的 精 
度 为 止 。 经 过 SPA 选择 特征 波长 后 ,入 选 的 建 模 的 
变量 有 限 , 且 在 此 基础 上 在 进行 PLSR 回 归 ,经 过 


10, 
SPA-MLR 预测 的 散 点 图 
8 3 
e | 
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R?=0.980 48 
2 RPD=10.501 3 
0 1 1 1 3 
0 2 4 6 8 10 
实测 值 /% 


的 土壤 盐分 预测 模型 


Fig.2 The prediction model of soil salt content based on VN 


表 6 


不 同 建 模 算法 预测 精度 比较 


Tab.6 Comparison of the prediction based on the different model 


SMLR 预测 SPA-MLR 预测 PLSR 预测 SPA-PLSR 预测 
数据 变换 2 2 2 2 
RMSEP Ral RMSEP Ra RMSEP Ra RMSEP Ra 
原始 光谱 0.388 6 0.953 2 0.311 8 0.970 7 0.556 9 0.916 0 0.311 8 0.970 7 
归 一 化 0.287 6 0.979 2 0.278 3 0.980 5 0.322 8 0.973 7 0.449 2 0.950 5 
移动 平滑 0.362 0 0.956 2 0.361 6 0.956 8 0.365 7 0.954 3 1.0515 0.589 4 
多 元 散射 校正 0.331 2 0.969 0 0.286 3 0.977 0 0.3447 0.9667 0.355 6 0.963 9 
变量 标准 化 0.331 4 0.969 4 0.304 9 0.974 1 0.343 5 0. 967 0 0.370 9 0. 960 8 
卷 积 平滑 0.388 8 0.951 5 0.307 9 0.968 9 0.615 7 0.884 5 / / 
一 阶 导 数 0.425 6 0.948 9 0.267 5 0.980 5 0.337 8 0.969 3 0.324 2 0.970 1 
二 阶 导数 0.435 6 0.949 4 0.341 6 0.969 6 0.400 4 0.958 7 0. 446 8 0.947 7 


0.05 的 水 平 检验 后 ,丢失 了 部 分 有 效 信 息 , 人 选 的 
主 成 分 个 数 也 随 之 减少 ,导致 模型 的 预测 精度 略 低 。 

就 PLSR 与 SMLR 预测 结果 比较 而 言 ,除了 导 
数 预 处 理 的 情况 ,SMLR 建 模 精度 要 优 于 PLSR 模 
型 ,这 可 能 是 由 于 逐步 回归 建 模 仅 有 少量 波段 的 光 
谱 吸 收 率 能 进入 模型 ,而 且 逐 步 回 归 也 可 以 减少 变 
量 间 共 线性 问题 ,模型 对 于 所 选取 的 特征 波长 信息 
利用 较 多 ,逐步 回归 所 选择 的 变量 更 能 代表 所 有 的 
波段 。 当 然 PLSR 也 能 够 处 理 自 变 量 较 多 的 情况 且 
能 降低 变量 之 间 的 共 线 性 ,但 在 本 文 的 实验 结果 中 
不 及 SMLR 建 模 效 果 好 。 因 此 ,SMLR 相对 于 PLSR 
可 以 更 好 地 构建 光谱 与 土壤 盐分 含量 之 间 的 关系 模 
AY 


4 结论 


利用 高 光谱 成 像 技术 获取 南 狂 典型 沙壤土 光谱 
信息 ,并 同步 测定 土壤 电导 率 , 通 过 SMR 和 SPA 算 
法 提取 特征 波段 ,结合 MLR 和 PLSR 建立 回归 模 
型 ,实现 了 对 土壤 电导 率 的 快速 检测 。 在 文中 所 列 
的 7 种 预 处 理 方 法 中 , 较 优 的 预 处 理 方法 是 归 一 化 ， 
多 元 散射 ,变量 标准 化 和 一 阶 导数 ;预测 模型 SPA- 
MLR 建 模 的 预测 精度 要 优 于 SMLR ,直接 使 用 PLSR 
建 模 的 预测 精度 要 优 于 SPA-PLSR , SMLR 和 SPA- 
MLR 模型 的 预测 效果 相对 PLSR 和 SPA-PLSR 较 理 
想 ,表明 相对 于 PLSR 方法 而 言 ,MLR 更 适合 建立 典 
型 沙壤土 电导 率 预测 模型 。 

本 文通 过 地 面 实测 土壤 光谱 和 电导 率 数据 , 综 
合 分 析 对 比 了 不 同方 法 在 新 疆 南 疆 阿 拉 尔 地 区 沙 壤 
土 盐分 检测 时 的 精度 和 效率 ,建立 的 SPA-MLR 预测 
模型 的 Ra = 0. 980 5, RPD = 11. 50, RMSEP = 
0.278 3 ,为 快速 7A ff He or WU PA Se A Be] E A FK 


提供 了 一 定 的 理论 参考 依据 。 此 外 ,新 疆 南部 地 区 
土壤 盐 渍 化 严重 ,本 研究 在 地 面 尺度 准确 获取 的 典 
型 沙壤土 的 盐分 敏感 波长 和 监测 模型 ,为 后 续 使 用 
卫星 和 无 人 机 遥感 进行 大 范围 的 有 效 监测 提供 了 理 
论 依据 和 实用 模型 ,为 构建 多 尺度 、 多 类 型 土壤 的 盐 
渍 化 遥感 监测 提供 基础 数据 和 技术 支持 。 但 是 不 同 
类 型 土壤 的 含水 量 .质地 结构 以 及 有 机 质 含量 等 理 
化 指标 不 同 ,光谱 特征 和 建 模 方法 也 可 能 存在 一 定 
的 差异 ,针对 不 同类 型 土壤 的 光谱 预 处 理 方法 .特征 
波长 选择 方法 和 建 模 方法 可 在 未 来 工作 中 加 强 对 比 
研究 。 


参考 文献 ( References ) 


[1] PANG G,TAO W,JIE L,et al. Quantitative model based on field- 
derived spectral characteristics to estimate soil salinity in Minqin 
County ,China[ J]. Soil Science Society of America Journal ,2014 ， 
78(2) :546 - 555. 

[2] ZHAO Z,TASHPOLAT T,ZHANG F,et al. Spectral characteristics 
of soil salt content in typical oasis of Tarim River’s middle reaches 
[J]. Journal of Natural Disasters ,2012 ,21 (5 ) :72 -78. 

[3] 吴 明 珠 , 李 小 梅 , 沙 晋 明 . 亚热带 红壤 全 氮 的 高 光谱 响应 和 反 
演 特征 研究 [J. 光谱 学 与 光谱 分 析 ,2013 ,33 (11) :3111 - 
3115. [WU Mingzhu ,LI Xiaomei,SHA Jinming. Spectral inversion 


models for prediction of red soil total nitrogen content in subtropi- 


cal region (Fuzhou) [ J]. Spectroscopy & Spectral Analysis ,2013 , 
33(11) :3111 -3115. ] 

[4] NAGY A,PETERRICZU,GALYA B,et al. Spectral estimation of 
soil water content in visible and near infra-red range[ J]. 2014 ,3 
(3) :163 -171. 

[5] ZHANG P,LI Y. Study on the comparisons of the establishment of 
two mathematical modeling methods for soil organic matter content 
based on spectral reflectance[ J]. Spectroscopy & Spectral Analy- 
sis 2016 ,36(3 ) :903 -910. 

[6] 张 晓 光 , 黄 标 , 季 峻 峰 , 等 . AET HY ULL Sb ies EE AY AR AG ER os 
土 盐分 定量 模型 研究 [J]. 光谱 学 与 光谱 分 析 ,2012 ,32(8): 
2075 - 2079. [ ZHANG Xiaoguang, HUANG Biao JI Junfeng, et 


al. Quantitative prediction of soil salinity content with visible-near 


201811.00063v1 


chinaXiv 


[10 


[11 


[12 


[14 


[15 


] 


] 


] 


wo 


Lo 


ChinaXiv 合 作 期 刊 


王涛 等 :MLR 和 PLSR 的 沙壤土 盐分 含量 光谱 检测 对 比 研 究 


infrared hyper-spectra in northeast China [ J]. Spectroscopy & 
Spectral Analysis ,2012 ,32(8) :2075 - 2079. | 

HRR, TEN, FHA, T. 基于 表 观 电导 率 与 实测 光谱 的 干 
早 区 湿地 土壤 盐分 监测 []]. 中 国 沙漠 ,2016,36(5):1365 - 
1373. [ YANG Aixia, DING Jianli, LI Yanhong, et al. Apparent e- 


lectronic conductivity and measured spectral for monitoring soil salt 
content in Arid Lakeside Wetland[ J ] . Journal of Desert Research, 
2016 ,36(5) :1365 - 1373. | 

GOLDSHLEGER N ,CHUDNOVSKY A, BEN-BINYAMIN R. Pre- 
dicting salinity in tomato using soil reflectance spectra[ J |. Interna- 
tional Journal of Remote Sensing ,2013 ,34(17) :6079 - 6093. 
NETO O R, TEIXEIRA A,LEAO R, et al. Hyperspectral remote 
sensing for detecting soil salinization using ProSpecTIR - VS aerial 
imagery and sensor simulation[ J ]. Remote Sensing ,2017, 9(1) :1 
- 16. 

WAS, EZR, ELIE, SE. 土壤 含 盐 量 与 电导 率 的 高 光谱 反 演 
精度 对 比 研究 [J]. 光谱 学 与 光谱 分 析 ,2014,34(2):510 - 
514. [ PENG Jie, WANG Jiaqiang, XIANG Hongying, et al. Com- 


parative study on hyperspectral inversion accuracy of soil salt con- 


tent and electrical conductivity | J]. Spectroscopy and Spectral A- 
nalysis ,2014 ,34(2) :510 -514. | 

届 永 华 , 段 小 亮 ,高 鸿 永 , 等 . 内 蒙古 河套 灌区 土壤 盐分 光谱 定 
st Od OT OF FE I]. 光谱 学 与 光谱 分 析 ,2009 ,29(5 ) :1362 - 
1366. [ QU Yonghua, DUAN Xiaoliang, GAO Hongyong, et al. 


Quantitative retrieval of soil salinity using hyperspectral data in the 


region of Inner Mongolia hetao irrigation district [ J ]. Spectroscopy 
& Spectral Analysis ,2009 ,29 (5 ) :1362 - 1366. | 

丁 建 丽 , 伍 漫 春 , 刘 海 霞 ,等 . 基于 综合 高 光谱 指数 的 土壤 
盐 渍 化 监测 研究 [ 相 . 光谱 学 与 光谱 分 析 ,2012 ,32(7) :1918 - 
1922. [DING Jianli ,WU Manchun ,LIU Haixia,et al. Study on the 


> 
= 


soil salinization monitoring based on synthetical hyperspectral index 
[J]. Spectroscopy & Spectral Analysis, 2012 , 32 (7): 1918 — 
1922. ] 

段 鹏 程 ,能 黑 钢 ,ZESR SR ,等 . 不 同和 干扰 程度 的 盐 渍 土 与 其 光谱 
反射 特征 定量 分 析 [ J]. 光谱 学 与 光谱 分 析 ,2017 ,37(2) :571 
—576. [ DUAN Pengcheng, XIONG Heigang, LI Rongrong ,et al. A 


auantitative analysis of the reflectance of the saline soil under dif- 


ferent disturbance extent | J]. Spectroscopy & Spectral Analysis, 
2017 ,37(2) :571 -576. | 

李 亚 莉 , 乔 江 飞 , 董 天 宇 ,等 .不 同 质地 盐 渍 化 土壤 水 盐 含 量 的 
uC IK [J]. 应 用 生态 学 报 ,2016,27(12 ) :3807 - 3815. 
[LI Yali, QIAO Jiangfei, DONG Tianyu, et al. Hyperspectral inver- 


sion of soil water and salt content in soils with different textures 
[J]. Chinese Journal of Applied Ecology ,2016 ,27 ( 12 ) :3807 — 
3815. | 

刘 娅 , 潘 贤 章 , 王 昌 昆 ,等 . 基于 可 见 - 近 红 外 光谱 的 滨海 盐 土 
土壤 盐分 预测 方法 [J]. 土壤 学 报 ,2012,49 (4):824 - 829. 
[LIU Ya,PAN Xianzhang, WANG Changkun, et al. Prediction of 


coastal saline soil salinity based on vis-nir reflectance spectroscopy 


[16] 


[17 


wo 


[18] 


[19] 


[20] 


[21] 


[22] 


[24] 


[J]. Acta Pedologica Sinica,2012 ,49 (4) :824 - 829. ] 

ZHAN X Y,ZHAO N,LIN Z Z,et al. Effect of algorithms for cali- 
bration set selection on quantitatively determining asiaticoside con- 
tent in centella total glucosides by near infrared spectroscopy[ J ] . 
Spectroscopy & Spectral Analysis ,2014 ,34( 12) :3267 - 3272. 
尼 珍 , 胡 昌 勤 , 冯 芳 . 近 红 外 光谱 分 析 中 光谱 预 处 理 方法 的 作 
及 其 发 展 []]. 药物 分 析 杂 志 ,2008 ,28 (5 ) :824 - 829. [NI 
Zhen,HU Changqin, FENG Fang. Progress and effect of spectral 
data pretreatment in NIR analytical technique[ J]. Chinese Journal 
of Pharmaceutical Analysis ,2008 ,28 (5) :824 - 829. | 

夏 俊 芳 , 李 培 武 , 李 小 垦 ,等 .不 同 预 处 理 对 近 红 外 光谱 检测 脐 
HE VC 含量 的 影响 [J]. 农业 机 械 学 报 ,2007,38(6) :107 - 
111.[XIA Junfang, LI Peiwu, LI Xiaoyu, et al. Effect of different 


pretreatment method of nondestructive measure vitamin C content 


of umbilical orange with near-infrared spectroscopy [ J |. Transac- 
tions of the Chinese Society for Agricultural Machinery , 2007 , 38 
(6) :107 -111. ] 

OD, HE WERE, EA BE, S.A EB E TEREX i BE AY 
盐 渍 土 盐分 含量 预测 [ 站. 光谱 学 与 光谱 分 析 ,2013 ,33(10 ) : 
2771 -2776. [| LIU Ya,PAN Xianzhang, WANG Changkun ,et al. 


Predicting soil salinity based on spectral symmetry under wet soil 


condition[ J]. Spectroscopy & Spectral Analysis ,2013 ,33 (10 ) : 
2771 -2776. | 
陈 文 倩 , 丁 建 丽 , 谭 娇 ,等 . 干旱 区 绿洲 植被 高 光谱 与 浅 层 土壤 
含水 量 拟 合 人 研究 [ J/OL]. 农业 机 械 学 报 ,2017, (9):1 -16 
[CHEN Wenqian, DING Jianli, TAN Jiao, et al. Analysis vegeta- 


tion of hyperspectral reflectance and shallow soil water content in 


arid area[ J |. Transactions of the Chinese Society for Agricultural 
Machinery ,2017 , (9) :1 -16. | 

PHAZE , 塔 西 甫 拉 提 “' 特 依 拜 , 买 买 担 ' 沙 吾 提 ,等 . 荒漠 土壤 
有 机 质 含量 高 光谱 估算 模型 [Jj. 农业 工程 学 报 ,2014 , 30 
(16): 113 - 120. [ HOU Yanjun, TIYIP Tashpolat, SAWUT 


Mamat ,et al. Estimation model of desert soil organic matter content 


using hyperspectral data[ J]. Transactions of the Chinese Society of 
Agricultural Engineering ,2014 ,30(16) :113 - 120. ] 

彭 翔 , 胡 丹 , 曾 文 治 , 等 . 基于 EPO - PLS 回归 模型 的 盐 渍 化 土 
壤 含 水 率 高 光谱 反 演 [可 ]. 农业 工程 学 报 ,2016 ,32(11 ) :167 - 
173. [ PENG Xiang, HU Dan,ZENG Wenzhi, et al. Estimating soil 


moisture from hyperspectra in saline soil based on EPO-PLS regres- 


sion[ J]. Transactions of the Chinese Society of Agricultural Engi- 
neering ,2016 ,32(11) :167 -173. | 

ZENG W Z,HUANG J S, XU C, et al. Hyperspectral reflectance 
models for soil salt content by filtering methods and waveband se- 
lection[ J]. Ecological Chemistry & Engineering S ,2016 ,23(1): 
117 - 130. 

FARIFTEH J, MEER F V D, ATZBERGER C, et al. Quantitative 
analysis of salt-affected soil reflectance spectra; A comparison of 
two adaptive methods (PLSR and ANN) [J]. Remote Sensing of 
Environment ,2007 ,110(1) :59 -78. 


A comparison of the salt content in sandy soil between the MLR 
model and PLSR model 
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Abstract: In order to monitor typical soil salt content (sandy loam soil) in South Xinjiang, China quickly and ef- 
fectively ,and to improve the precision of the soil salt content estimation model by removing the noise of soil hyper- 
spectral absorbance , this paper focused on the inversion relationship between soil spectrum and electrical conductivi- 
ty (EC) by using multiple spectral pretreatment methods, and then the multiple linear regression (MLR) and the 
partial least squares regression (PLSR) modelling were applied to establish the salt content model based on the hy- 
perspectral analysis technique. The effective and predictive capacities of different models were validated. This study 
took the typical arid area in South Xinjiang as the research object, obtained the hyperspectral data and EC by using 
Near-infrared spectrometer ( Zolix Gaia Sorter) and conductivity meter (DDS - 307) ,142 soil samples at 0 ~ 20 cm 
depth were collected and these samples were highly representative for the EC. Seven pretreatment methods were 
used to pretreat the original spectral data, such as vector normalization ( VN ) , multiplicative scatter correction 
(MSC) , standard normal variate ( SNV ) , moving-average (MA) smoothing, Savitzky-Golay (SG) smoothing, first 
derivative (1 - Der) and second derivative (2 - Der) ,then the characteristic wavelengths were extracted with step- 
wise multiple regression (SMR) and successive projection algorithm (SPA) , which were used as input variables of 
MLR and PLSR modeling. The results showed that the optimum pretreatment methods were VN,MSC,SNV and 1 — 


2 


Der. According to different pretreatments ,in the stepwise multiple linear regression (SMLR) prediction model R, 
was greater than 0.95 ,RPD was greater than 6. 2 ,and RMSEP was less than 0. 44. In the multiple linear regression 
model based on the successive projections algorithm (SPA-MLR) R,,,” was greater than 0.96 ,RPD was greater than 
6.9,and RMSEP was less than 0.36, which were better than those in SMLR. In the partial least squares regression 


prediction model R, was greater than 0. 88 ,RPD was greater than 4.4,RMSEP was less than 0. 62 and in the par- 


2 


tial least squares regression model based on the successive projections algorithm (SPA-PLSR) R,, was greater than 


0.59,RPD was greater than 2.4,and RMSEP was less than 1. 1 which were less than those in SMLR and SPA- 
MLR. The best prediction of SMLR and SPA-MLR after vector normalization (VN) was as the follows: RMSEP =0. 


287 6,R =0. 979 2,RPD =9. 907 8 and RMSEP =0. 278 3,R,,, =0. 980 5,RPD =11.50 which had fewer 


val val 


characteristic wavelengths. So the VN is a best effective pretreatment method. It is more suitable to establish the pre- 
diction model of soil conductivity in jujube tree by MLR with PLSR. It could be an important part in future resear- 
ches how to choose the right algorithm to analyze the soil salt content with the spectral data. 

Key words: soil electrical conductivity; multiple linear regression; Partial least squares regression ; 


Hyperspectral 


